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Where have all the black holes gone? 
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Abstract. We have calculated stationary models for accretion disks around super-massive black holes in galactic nuclei. Our 
models show that below a critical mass flow rate of ~ 3 10 -3 MEdd advection will dominate the energy budget while above 
that rate all the viscously liberated energy is radiated. The radiation efficiency declines steeply below that critical rate. This 
leads to a clear dichotomy between AGN and normal galaxies which is not so much given by differences in the mass flow rate 
but by the radiation efficiency. At very low mass accretion rates below 5 10 JI4dd synchrotron emission and Bremsstrahlung 
dominate the SED, while above 2 10 _4 MEdd the inverse Compton radiation from synchrotron seed photons produce flat to 
inverted SEDs from the radio to X-rays. Finally we discuss the implications of these findings for AGN duty cycles and the 
long-term AGN evolution. 
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1 . Evidence for Black Holes and 
Advection-Dominated Accretion 

The existence of black holes (BH) in the centers of galaxies 
is now widely accepted and the best mass determinations are 
known for Sgr A* in the Galactic Center (M B h = 2.6 10 6 M Q ) 
from the stellar velocity dispersion, and for NGC4258 from 
keplerian rotation of maser spots in an accretion disk. Beside 
these very low luminosity AGNs, the masses of BH in quasars 
have been estimated from their continuum luminosity and the 
H/3 line width (Laor 2000 ). It turns out, that radio-loud quasars 
and radio galaxies, which show powerful jets, are found in large 
elliptical galaxies with the most massive BHs ~ 10 9 A/ Q . 

Seyfert Galaxies, on the other hand, are spiral galaxies 
which host an AGN . Activity in the nucleus, which is pow- 
ered by accretion into a BH, can be discriminated against star- 
bursts in Seyferts from radio and X-ray observations. Radio 
cores and jets with brightness temperatures ab ove 1 8 K have 
been detec ted in some S eyfer ts (Ulvestad et al. 1999 ; Mundell 
et al. |2000| ; Falcke et al. |200C| ). Their flux stability over several 
years exculdes radio supernovae as the power supply. The X- 
ray emission shows rapid variability and in some cases a red- 
shifted Fe Ka line, which is an indicator of relativistic mo- 
tion in the accretion disk around the BH. The masses of BHs 



Send offprint requests to: T. Beckert, Max-Planck-Institut fiir 



Radioastronomie, Auf dem 
tbeckert@mpifr-bonn.mge.de 
Correspondence to: W.J. 
Astrophysik, Tiergartenstr. 
wjd@ita.uni-heidelberg.de 



Hiigel 69, 53121 Bonn, Germany, 

Duschi, Institut fiir Theoretische 
15, 69121 Heidelberg, Germany; 



in some Seyfert 1 galaxies have been measured by reverber- 
ation mapping of variable an d corr elated continuum and line 
emission (Peterson & Wandel 200C ). These measurements are 
in reasonable agreement with the Mbh-c relation of enclosed 
mass (~ Mbh) versus velocit y disp ersions a in bulges of nor- 
mal galaxies (Gebhardt et al. |2000| ). It is therefore reasonable 
to assume the existence of supermassive BH in most elliptical 
galaxies and spirals with bulges. 

While in the high luminosity objects both jet and accre- 
tion disk can be identified in the spectrum, the situation is dif- 
ferent in less luminous AGNs like weak Seyfert Galaxies and 
LINERs. But even here small scale jets are commonly found 
(Falcke et al. |2000| > and argue for the existence of BHs. For 
instance NGC4258 is an interesting transition object show- 
ing both an outer irradiated thin accretion disk and a small 
scale radio jet. A geometrically thin standard accretion disk 
close to the black hole can not be identified, but the ioniz- 
ing X-rays maybe produced at the base of the jet, which can 
be identical with the proposed advection-dominated accretion 
flow (ADAF) within 100i?s (Gammie, Narayan & Blandford 
1999; Rs'. Schwarzschild radius). At even lower luminosities 
the Galactic Center (Sgr A*) with a BH mass of 2.6 1O 6 M 



(Genzel et al. 1997) is the only visible AGN with a power 
output of w 2 10~ 10 LEdd- A comparable object in any other 
galaxy (spiral or elliptical) would not be seen as an AGN. 
Assuming an spherical and adiabatic Bondi inflow, cooled 
only by Bremsstrahlung, gives an radiation efficiency r\ = 
9 10- 3 (i/LEdd) (Frank, King & Raine|1992[). Bremsstrahlung 



will be emitted in X-rays and the Chandra detection (Baganoff 
et al. 2001a) of Lx ~ 2 10 33 erg s _1 is consistent with a mass 
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Fig. 1. Spectral energy distribution (SED) for 
ADAF models with a black hole mass of 
2.610 6 Ajf Q . The mass accretion rates are M = 
(1(T 6 , 2 1CT 6 , 5 1CT 6 , 1(T 5 , 2 1CT 5 , 5 1(T 5 , 10~ 4 , 1.27 1(T 4 ) 
yr _1 . The scale-free accretion rate for this BH is 
m — 17.33[M/(M yr _1 )]. The spectral energy flux in- 
creases for increasing mass accretion rate with an exception at 
M = 1.27 1CP 4 , which is plotted with a dashed curve. 



accretion rate of M = 1.6 10 _6 Mq yr _1 and an efficiency of 
r\ = 2.2 10~ 8 . The sub-mm luminosity of Sgr A* is about 30 
times larger than the X-ray flux and makes Sgr A* a unique ob- 
ject. We will discuss a specific ADAF model for Sgr A* in Sec. 
0. The Bondi flow faces at least two problems: it does not allow 
for any possible angular momentum of the inflow and does not 
include magnetic fields, which lead to synchrotron emission at 
radio frequencies and synchrotron self-compton cooling. Both 
can be accounted for in ADAF models. They provide a reason- 
able explanation for the spectral energy distribution (SED) of 
Sgr A* with a mass accretion of « 1O _6 M0 yr _1 and a radia- 
tive efficiency of 1CP 5 . Beside the basically unresolved radio 
core of Sgr A*, it is not possible to identify a jet in the Galactic 
Center. 

In this paper we will explore the hypothesis, that most of 
the normal galaxies without substantial AGN activity contain 
supermassive black holes M > 10 6 Af© some of which have 
been active during the quasar phase 0.3 < z < 5 (z being the 
cosmological redshift) and are quietly accreting in an ADAF 
mode today. Spectral properties of ADAFs with rather large 
mass accretion rates are explored in Sec. 2. The total luminosi- 
ties and sp ectral energy distributions are of interest for weak 
AGNs (Ho 1999 ). We investigate the transition from standard 
thin disk accretion to ADAFs and vice versa as an upper limit 
in the mass accretion rate for ADAFs in Sec. fjj. The combined 
consequences for accretion in normal galaxies are discussed in 
Sec.|. 
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Fig. 2. SED for an ADAF around a 10 9 Af o black hole. 
The mass accretion rate is m = 3.6 10~ 4 corresponding to 
8 lO _3 M0/yr. The seed photon flux is represented as dashed 
line, the first and second Compton humps are dash-dotted and 
dotted respectively. Further Comptonization is important but 
^individual humps are smeared out. Only the total SED is shown 
as solid line. The dependence of the spectral shape on the black 
hole mass is weak, while the total flux scales with the mass ac- 
cretion rate according to Eq. (0). 



2. Radiation characteristics of ADAFs 

The calculation of emission spectra from ADAFs are based 
on self-similar solutions for density, inflow velocity, rotation, 
and ion temperature presented in Beckert (2000), which are 
a generalisation of the Narayan & Yi (1994) solution in the 
Newtonian limit. Gravitational redshift and the small volume 
of the general-relativistic region makes it reasonable to extend 
the solution down to the black hole horizon at i?s- Aspects 
of the Kerr metrics are not considered. The electron tempera- 
ture is calculated from thermal balance between Coulomb heat- 
ing by ions, adiabatic compression, advective energy transport, 
and radiation cooling. The primary radiation mechanisms are 
Bremsstrahlung and synchrotron emission. Both are local in 
the sense, that they only depend on the local state of the flow. 
At large mass accretion rates multiple inverse Compton scatter- 
ing (IC) of synchrotron photons becomes the dominant cooling 
mechanism. The photon flux to be scattered is produced at dif- 
ferent radii in the flow and is highly anisotropic, which makes 
it a non-local process. Inverse Compton scattering and the re- 
sulting cooling rate modify the density and temperature of the 
ADAF and are treated iteratively. For calculating the IC radi- 
ation, we used a method proposed by Poutanen & Svensson 
(1990 for thermal electrons, which we have modified to al- 



low for the quasi-spherical geometry of ADAFs. From the total 
cooling we take the global radiation efficiency £R a d = Q~/Q + 
to recalculate the self-similar ADAF structure and electron 
temperature. We then seek convergence of the assumed and 
posteriori calculated £R a d to get globally consistent ADAF 
models with correct radiation spectra. For the ADAF we as- 
sume energy equipartition between ions and magnetic field. 
This determines the synchrotron emission and provides addi- 
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tional pressure to support the flow and lower the adiabatic in- 
dex below 5/3 to make ADAFs possible at all. For the vis- 
cosity we use the standard a-prescription with a = 0.1 well 
above the critical value of a cr it ~ 1CP 2 for the transition to 
convection-dominated accretion flows (CDAFs) described by 
Narayan, Igumenshchev & Abramowicz ( 2000) ). No outflows, 
jets or wind infall as investigated in Beckert (200C) are as- 
sumed, even though a wind infall is suggested for the Galactic 
Center (Melia & Coker 1999). The radial viscous break due to 



bulk and shear viscosity is included, but it does not dominate 
the dynamical state of the flow for a < 0.1. 

From the described model we can construct spectral en- 
ergy distributions (SED) for different mass accretion rates. The 
mass of the central BH has only a weak influence on the SED, 
which is not included in the scaling of M to the Eddington ac- 
cretion rate, and is not considered here. We assume a mass of 
2.6 10 6 M Q for the SED in Fig.[l] appropriate for the Galactic 
Center, with a corresponding Eddington limit of Mecm = 
0.0577M© yr" 1 (in this paper we define the Eddington accre- 
tion rate with an efficiency r\ = LEdd/(AfEddC 2 ) = 0.1) and 
use the scale free accretion rate rh = Af/AfEdd in the fol- 
lowing discussion. The spectral luminosity in Fig. [j] scales as 
the black hole mass, vL v oc Mbh- For comparision we show 
in Fig. | the SED for an ADAF arround a 10 9 M o black hole 
with to = 3.6 10 -4 . Fig. ^Jdemonsttrates that the radiation ef- 
ficiency for this flow is independend of the black hole mass. 

The SEDs for rh between 1.7 10~ 5 and 2.2 10~ 3 are shown 
in Fig. |TJ The presented model spectra are accurate above 30 
GHz, and they show that the synchrotron emission rises in flux 
from 10 33 to 10 37 erg s _1 and shifts in frequency from 2 10 11 
Hz to 2 10 12 Hz at rh = 3.5 10~ 4 and back to smaller fre- 
quency for larger to. Above rh — 1.5 10~ 4 the Thomson opti- 
cal depth for synchrotron photons from central regions around 
3i?s is significant, and Compton scattering broadens the syn- 
chrotron peak and make it less prominent, compared to the 
IC emission. The dominating peak in the IC part of the spec- 
trum, which can be identified between the synchrotron and the 
Bremsstrahlung peaks at 10 19 20 Hz, is the second Compton 
peak of twice scattered synchrotron photons. The synchrotron 
seed photons are produced in a region closer to the BH than the 
Compton scattered radiation. So the seed photon flux for first 
Compton scattering is anisotropic, and most synchrotron pho- 
tons are scattered back into a high density and high tempera- 
ture region with the largest optical depth. The second Compton 
peak is therefore the dominant one, and the asymmetry between 
even and odd scattering order decreases thereafter, because the 
photon field to be scattered becomes more and more isotropic. 
The SED becomes flat or inverted due to multiple IC scatter- 
ing above 10 14 Hz for rh > 3 10~ 4 . The Bremsstrahlung peak 
will only be recognized below rh < 1.5 10 -4 . The peak does 
not shift very much in frequency as the maximum of electron 
temperature only varies between 2.5 10 9 K and 8 10 9 K, where 
the highest T e are achieved at to « 2 10 -4 very close to the 
horizon and the photons from that radius are significantly red- 
shifted. 

One major prediction of ADAF models is the different evo- 
lution of observable flux in different frequency bands. This is 
expressed in the scale-free spectral luminosity l v . We scale the 
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Fig. 3. Normalized flux l v (see Eq. |Tj) at 86 GHz (solid line), in 
K (dash-dotted) and V band (dotted), at 1 keV (dashed) and 100 
keV (3 x dotted-dashed) as a function of mass accretion rate. 
The models are calculated for a 2.6 10 6 M Q black hole. The 
efficiency r\ in Eq. ([!]) is 10% and the SED is assumed to peak 
atlOOkeV. 

radiation flux L v to the total, frequency integrated luminosity 
of the specific model from Eq. (0) and define 



r)M Edd c 2 



m 



10" 



-2.3 



llVr, 



100 keV 



(1) 



We see in Fig.^J that the importance of synchrotron emission 
at 86 GHz decreases with rising to, the X-rays follow the to- 
tal luminosity at first, and become more important when they 
are dominated by IC emission above rh > 2 10~ 4 . The im- 
portance of IC emission is most prominently seen in K and 
V band. The flux is rapidly rising with to and saturates above 
to = 2 10~ 4 when IC is dominating the cooling and therefore 
IC follows the total luminosity. The described spectra show 
only a weak dependence on the absolute mass of the BH be- 
tween 10 4 and 10 9 M Q , for which we have tested the mod- 
els. Larger BH masses imply smaller densities in the accre- 
tion flow at the same to. Consequently the electron temperature 
also becomes smaller due to weaker Coulomb coupling. As the 
only significant consequence for the SED, the position of the 
synchrotron peak is anticorrelated with Mbh, and it shifts to 
smaller frequencies for larger BH masses. 

3. A Note on the Galactic Center Source Sgr A* 

The enigmatic radio source Sgr A* in the Galactic Center is 
coincident with the center of gravity of an enclosed mass of 
2.6 10 6 Mm. It was considered to be one of the test cases for 



ADAF models (Narayan et al. 1998 ), but the SED of the source 



poses three problems for standard ADAFs. (1) The observed 
radio spectrum is much flatter than predicted, so that only the 



sub-mm bump (Falcke 1999) is nowadays attributed to the ac- 



cretion flow. Most of the radio emission at cm-wavelength must 
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Fig. 4. Spectral energy distribution (SED) for an ADAF model 
with a wind infall appropriate for Sgr A* . The black hole mass 
is 2.6 10 6 Mq and the mass accretion rate at the last stable orbit 
M = 9 1O -7 M0 yr" 1 corresponding to ra = 1.6 10~ 5 . This 
special model is briefly described in Sec. ||. The dashed line 
shows the synchrotron and bremsstrahlung emission. The first 
three inverse Compton stages are given as dash-dotted lines. 
The total spectrum (thick line) has a remarkably flat radio slope 
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consistent with the observations. 



then be produced by an outflow or jet (Falcke & Markoff 2000). 
(2) The X-ray s pectrum as derived from Chandra observations 
(Baganoff et al. 2001b ) has a different slope than expected from 
thermal bremsstrahlung coming from the ADAF. ( 3) The ob- 
served rapid variability in X-rays (Baganoff et al. 2001a ) re- 
stricts the size of the variable emitting region to less than 10i?s- 
The spectrum at high X-ray fluxes is harder than at low flux 
levels. This can be explained by inverse C ompto n emission of 
relativistic electrons in a jet (Markoff et al. 2001 ) but even a jet 
has to be powered by an accretion process and bremsstrahlung 
emission of the accreting gas is unavoidable. In contrast to 
these recent scenarios, here we present an ADAF-wind infall 
model, where the gas in the accretion flow is heated by wind 
infall at all radii (Beckert, 2000) with a steeper density pro 



file Ear 1/2 13 than normal ADAFs. The synchrotron emis- 
sion dominates the SED (Fig. due to a strongly magnetised 

ADAF (3 P = PGas/Ptotal < 0.5. 

The wind infall is assumed to be strong (3 — 0.24 and ro- 
tates with ft/2 of the ADAF. The flow is magnetically domi- 
nated with (3p = 0.35 and the radiative efficiency e = 9 10 -5 
is larger than in the other models presented in this paper 
due to the larger electron temperature and the stronger mag- 
netic fields, which leads to increased synchrotron emission. 
Viscosity is described by an a-parametrisation with a — 0.08, 
lower than for the other ADAF models in Sec. 2, but convection 
is still expected to be unimportant (Narayan, Igumenshchev & 
Abramowicz 2000| ). This model gives a good fit to the radio 
spectrum but the problem with the X-ray observation persists. 



4. Limits on the mass accretion rate 

The luminosity of the synchrotron peak relative to 
Bremsstrahlung and IC emission depends only on m, 
and it is unaffected by the absolute mass of the BH. The global 
radiative efficiency £R a( j, defined in the previous section, gives 
the cooling rate Q~ with respect to the viscous heating rate. 
The heating rate itself depends on £R a( j in the ADAF model, 
and it is therefore not useful for comparing different models. 
In describing the efficiency of an accretion flow, it is more rea- 
sonable to use e = Q~ / \rjMc 2 ). In the following we assume 
a standard efficiency of r] = 0.1 . The quantity e gives the ratio 
between the actual luminosity of the accretion flow and the 
value expected for a thin and effectively radiating accretion 
disk at the same mass accretion rate. The radiative efficiency e 
is shown for our model calculations in Fig. [|. As described in 
the previous section, different emission mechanisms dominate 
at different accretion rates, and we do not expect e to follow a 
simple power law in m. The fit 



(340 m) 



2.2 



(2) 



therefore does not represent a theoretically derived physical 
law, and has to be taken with caution when used outside the 
scope of the calculated models. Nonetheless the results allow 
us to conclude that no consistent ADAF models with mass ac- 
cretion rates larger than m cr i t = 2.95 10 -3 are feasible. The 
reason is that IC emission increases rapidly with electron tem- 
perature and density. For increasing mass accretion rates both 
density and surface density increase linearly as long as the ac- 
cretion velocity stays the same (this is true for ADAFs with 
e < 0.3, which are consequently advection-dominated). The 
maximum T c in the flow depends only weakly on m for IC 
dominated cooling. The Thomson optical depth is 



(3) 



and the mean energy which a photon gains from one scatter- 
ing is 2"f 2 hv w Q.hhv. The seed photons have v « 10 12 Hz, 




and upscattering stops if (27 



2\n 



jm e c r/(hv), which re- 



sults in photons gaining energy in up to n « 10 consecutive 
scatterings. Starting from a synchrotron luminosity L s the IC 
luminosity will be L\ C ps (2^ 2 T) n L s , if 2j 2 t > 1, and the 
spectrum will be flat for 2j 2 t = 1. In the case of 77 = a = 0.1 
and r ~ 3i?s this criterion is fulfilled, if m = 1.3 10~ 3 . This 
is a reasonable order of magnitude estimate, when compared to 
the model spectra described above. With a constant synchrotron 
luminosity of 5 10 _7 LEdd> we ^ n< ^ m tn e IC dominated regime 



Lie = 510" 



3 ) n -^Edd 



An upper limit for the mass accretion rate is 

to c = 0.5377a, 



(4) 



(5) 



which is in rough agreement with the upper limit found by 
fitting the total luminosity in Eq. (||). It must be noted that 
this limit m c has a different a-dependence than other lim- 
its m cr it = £a 2 (Narayan 1996); Esin, McClintock, Narayan 
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Fig. 5. Radiation efficiency e as a function of mass accretion 
rate m with efficiency 77 = 0.1 to compare with the expec- 
tation for efficient radiation cooling. Model calculations are 
made for a 2.62 10 6 M Q black hole (symbol: stars), for 10 8 M Q 
(diamonds), and for a 1O 9 M (triangles) black hole. The es- 
timated errors for the numerical convergence of assumed and 
posteriori derived e are given for the 2.6 10 6 Af Q black hole. 
The data are fitted with Eq. fy. 



1997) which give a wide range for £ between 0.1 and 1.3 de- 



pending on details of the models. 

Another local criterion for the existence of ADAFs is given 
by the imbalance of Bremsstrahlung cooling and viscous heat- 
ing. At large radii the electron cooling rate, which is coupled to 
the ion heating rate by the radiation efficiency, decreases faster 
than Coulomb coupling between thermal ions and electrons. 
The electron temperature is therefore close to the ion temper- 
ature at these radii. The ions are close to the viral temperature 
as long as the radiative efficiency is significantly smaller than 
1, which is the case for all calculated models here. In the re- 
gion where T c and T lon are equal, the Bremsstrahlung cooling, 
which dominates at large radii, decreases as r~ 5 / 2 , while the 
viscous heating falls off as r -3 . At a critical outer radius, the 
radiative efficiency is 1 and no ADAF is possible at larger radii. 
This outer radius is found to be 



5. Discussion - Consequences for AGN 
Evolution 

We find that the radiation efficiency depends strongly on the 
mass accretion rate (oc to 2 - 3 ) below a certain critical value 
TO cr it of the mass accretion rate (in units of its Eddington 
value). Above m cr it, e = 1. In other words, below the criti- 
cal mass accretion rate, the efficiency decreases so rapidly that 
one can discern two regimes: 



s = 



1 for TO > TO cr ;t 
for to < m cr jt 



(7) 



with only a small transition zone in to below m cr it where e 
differs significantly from 0. For the purpose of our present dis- 
cussion, however, the approximation of Eq. (Q) suffices. 

Furthermore, our numerical models show that the relation 
outlined by Eq. (^) holds for the entire range of Mbh investi- 
gated. In the following we will use the value m crit — 0.003 as 
derived by interpolating our numerical models (Eq. fy. 

5.1. The onset of nuclear activity 

The existence of m cr it with the properties discussed above 
translates into a fairly sharp transition between an active and an 
inactive state of a galaxy. As soon as to falls below m CIit the 
radiation efficiency of the accretion decreases dramatically. In 
other words, already a relatively small change in the mass flow 
rate around TO cr ;t suffices to "switch off an AGN, and vice 
versa. The difference between a normal and an active galaxy is 
then not due to a difference in to which is as large as the dif- 
ference in luminosities between the two classes. A much more 
important reason is the steep decline in the radiation efficiency 
for the accretion rates below which the disk turns advection- 
dominated. 

Our numerical models predict m cr it ~ 0.003. This is in 
good agreement with observations (e.g., Peterson & Wandel 
[200C , who find AGN only in the luminosity range between 
~ 10~ 3 and 1 of the Eddington luminosity, or — equivalently — 
the Eddington mass accretion rate). In our interpretation the 
lack of galactic nuclei below 10~ 3 ^Edd is not due to a lack of 
galaxies with mass accretion rates below 10~ 3 AfEdd but rather 
due to the steep decline of radiation efficiency below this criti- 
cal value. 



R out = 5 10V m- 1 R s , (6) 

which is valid for a < 0.1 and means that for to > 0.13 
Bremsstrahlung from an ADAF is more efficient than the vis- 
cous heating outside the marginally stable orbit at R = 3Rs- 
From studies of the Galactic Center we know that no stan- 
dard thin accretion disk exists within 10 5 i?s an d the mass 
accretion rate from spectral fitting of an ADAF model gives 
to w 1.6 10~ 5 . From Eq. (^) this results in a lower limit esti- 
mate of a > 0.01. 



5.2. The general properties of the evolution of AGN 
and of their black holes 

For a supermassive black hole of mass Mbh, one can give an 
average accretion rate M over its age tbh of 

M<^ (8) 
tbh 

We can only give an upper limit of M, because there is the pos- 
sibility of a non-negligible seed mass of the black hole which 
is not due to this type of accretion processes. 



6 



T. Beckert & W.J. Duschl: Where have all the black holes gone? 



If, in addition, we assume, that the age of the BH is not 
very much shorter than the age of its host galaxy and thus the 
Hubble time (at the location of the black hole), th, we can write 

M<^£ (9 ) 

TH 

In the following, we express the time in units of 10 10 yr, and the 
mass flow rates in units of the Eddington rate. We introduce the 
abbreviations rh = M /AfEdd, and no = r/10 10 yrs, and get 
an average mass accretion rate of 

m < 0.0477— . (10) 

no 

In terms of the Eddington accretion rate, the average accretion 
rate rh is independent of the accreting black hole's mass. For a 
constant accretion efficiency 77, it is m oc a function of the 
age of the black hole, or — for that matter — the Universe. This 
leads to two interesting consequences: 

- In terms of the Eddington rate, the accretion rate declines 
as the black hole, the galaxy, and the Universe as a whole 
evolve; 

- The present-day average accretion rate is around 10~ 3 of 
its Eddington value. 

The present-day average accretion rate, m/Q) is fairly close 
to the critical accretion rate 7h cr it ~ 0.003 given by s = 1 
(Eq. ||). For 771 > m C rjt the radiation efficiency is unity, while 
for smaller rh it decreases sharply oc m 2 3 . This means that 
relatively small changes in the momentary accretion rate are 
capable of transferring a galaxy from an inactive to an active 
state and vice versa. In the course of the further evolution of 
the Universe, it will become harder and harder for galaxies to 
turn active as the average accretion rate (in Eddington units) 
becomes smaller and smaller. This is compounded with and 
strengthened by a diminishing supply of material available for 
accretion, i.e., by an additional decrease of the absolute value 
of the average accretion rate. Extrapolating back to earlier cos- 
mological epochs, we find (because of smaller tbh and higher 
771) that the likelihood for a galaxy to be active was higher on 
two grounds: (1) The larger rh the smaller a fluctuation suffices 
to turn a galaxy active. (2) The further we go back in the evolu- 
tion of the Universe, the larger was the supply of gas available 
for accretion. 

5.3. Constraints on the AGN duty cycle 

As discussed above, the mass Mbh and age tbh of a black 
hole define an upper limit for its average mass accretion rate. 
During phases of activity, the mass flow rate must be larger than 
?Ti. Let us — for the purpose of a crude estimate — assume that 
we have only two states, namely the AGN phase, characterized 
by a mass flow rate ttiagn > ihcxit lasting for a period of 
time of 6th, and a normal galaxy phase for which the mass 
flow rate 77i norma i < rh cr jt is correspondingly smaller so as to 



maintain the average value in. Let us, moreover, assume[| that 

(1 — 9)rh nolma i -C 0771 agn, then we get for the duty cycle 

777. — 77l norma l 771 

9 = > . (11) 

TOAGN ~ ra normal 2m A GN 

toagn > w cr it needs to be fulfilled for an activity phase to 
occur at all. The stronger the activity of a galaxy is the shorter 
it can be maintained. For instance, for an AGN operating at 
its Eddington limit, i.e., at 771agn = 1 this means that it can- 
not stay at this level of activity — integrated over all individual 
phases of activity — for longer than a fraction 9 of its entire evo- 
lution, i.e., some 10 7 yr in the present-day Universe. A super- 
Eddington activity level can be maintained only for an even 
shorter period of time. 

Derivation of a more detailed luminosity evolution of an 
AGN sample requires, however, a treatment more detailed than 
the above order-of-magnitude estimates. In particular it has to 
be investigated whether real-world galaxies can maintain a suf- 
ficiently high supply of matter for the accretion process over 
long enough a period of time. This then involves, for instance, 
questions about the accretion time scales and the long-term de- 
velopment of mass reservoirs. This topic, however, is beyond 
the scope of the present paper and will be addressed separately 
(Duschl & Strittmatter, in prep.) 

6. Conclusions 

We have shown that advection-dominated accretion flows into 
black holes display changing spectral energy distributions for 
different mass accretion rates. At very low mass accretion 
rates below M — 5 10~ 5 AfEdd synchrotron emission and 
Bremsstrahlung dominate the SED. Above M — 2 10 _4 MEdd 
the inverse Compton radiation from synchrotron seed pho- 
tons produce flat to inverted SEDs from the radio to X-ray 
bands. Multiple inverse Compton scattering is the most rele- 
vant cooling process for mildly relativistic electrons in two- 
temperature ADAFs, whenever the actual radiation efficiency 
in eq. (§) is larger than 0.07 or the total luminosity larger than 
4 10 _7 LEdd- The rapidly increasing cooling efficiency of the 
inverse Compton process sets an upper limit for the mass accre- 
tion rate of 77i cr it = 3 10~ 3 for ADAFs. At larger mass accre- 
tion rates radiative cooling balances the local viscous heating 
and advective energy transport is unimportant. The resulting 
flows are less hot, disk-like, and therefore much denser than 
ADAFs. Below the critical accretion rate, the radiation effi- 
ciency declines rapidly e oc 771 2 3 and the total luminosity is 

Ladaf = 3.45 10- 3 (77i/0.003) 3 - 3 i E dd • (12) 

In addition to the changing spectral behavior, the combi- 
nation of a transition value between ADAFs and normal ac- 
cretion disks (rhcrit ~ 3 10 -3 ), and the very steep decline of 
the ADAF radiation efficiency below this transition value to- 
wards smaller accretion rates (e oc rh 2 3 ) leads to an apparent 

1 This assumption is not necessary for activity to set in, as we 
have seen in the previous Sect. For the following, however, it is a 
handy assumption which does not influence the results of our order- 
of-magnitude estimates. 
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dichotomy in what one observes: Above m cr it> one will rec- 
ognize the sources as AGN; while below m cr it> practically no 
nuclear activity is observable. In other words, accretion into a 
super-massive black hole should be prominently visible only 
at mass flow rates above TO cr ; t , which is in good agreement 
with observed AGN distributions. The major difference be- 
tween AGN and normal galaxies is not so much a difference 
in mass flow rate — at least not by as much as the difference in 
luminosities may make us think — but rather one in the radiative 
efficiency. 

At the same time, this means that — at accretion rates 
around m C rit — small changes in the mass flow rate are suffi- 
cient to cause a strong difference in radiation efficiency and 
thus nuclear luminosity. In other words, the crossing of ra cr it 
acts almost like a switch which turns AGNs on and off. 

Finally, the combination of the black holes' masses and the 
mass accretion rate allowed us to put constraints on the duty 
cycle of AGN. It turned out that the most active AGN can main- 
tain this level of activity only for rather short time scales (of the 
order of some 10 7 years). 
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